atures is involved during the impulsive phase, and single tem- 
perature analysis is inadequate, given the improvements in 
the observations. Multi-thermal analysis techniques, 
however, are limited in the detail with which the tempera- 
ture distribution can be determined, even with more accurate 
observations. fjO-Jt _ X 9336 { 

5.3 THE ENERGETICS OF THE 
GRADUAL PHASE 


K.T. Strong, R.D. Bentley, P.L. Bommann, M.E. Bruner, 
P.J. Cargill, J.G. Doyle, J.R. Lemen, R. Pallavicini, 
G. Peres, S. Serio, G.M. Simnett, J. Sylwester, and N.J. 
Veck 


we compare our results with those in the chapter by Moore 
et al. (1980), who reached five main conclusions about the 
gradual phase: 

• the typical density of the soft X-ray emitting plasma 
is between 10 11 and 10 12 cm“ 3 for compact flares and 
between 10 10 and 10 11 cm -3 for a large-area flare; 

• cooling is by conduction and radiation in roughly equal 
proportions; 

• continual heating is needed in the decay phase of two- 
ribbon flares; 

• continual heating is probably not needed in compact 
events; 

• most of the soft-X-ray-emitting plasma results from 
‘ 4 chromospheric evaporation ’ ’ . 


5.3.1 Introduction 

The gradual phase of a flare is best characterized by the 
smooth rise and decay of the soft X-rays (see Figure 5A.2). 
Just from inspection of the soft X-ray light curves of some 
typical flares, it is evident that the gradual phase comprises 
several different stages. During the rise of the X-ray emis- 
sion, a completely different energy budget from that of the 
decay will be produced. Many flares show multiple peaks 
in their light curves, indicating that energy is released after 
the initial burst. The location of any such secondary energy 
release will also profoundly affect the relative importance 
of the terms in the energy budget of the flare. 

The soft X-ray emission produced by the higher- 
temperature ( > 10 7 K) plasma rises and decays more rapidly 
and peaks earlier than that produced by the cooler plasma 
(10 6 - 10 7 K). The impulsive phase usually occurs during 
the rise of the soft X-ray emission but there is no clear evi- 
dence of an impulsive component in the soft X-ray signal. 
Hence, in an energetics study the important question arises 
of what, if any, is the link between the impulsive and grad- 
ual phases. 

We have avoided using the term “thermal phase” be- 
cause it is misleading; a number of non-thermal energetics 
terms are associated with the gradual phase. The mass mo- 
tions are an example; it is not clear whether they are peculiar 
to the early stages of a flare or whether they are also present 
later in the flare, which seems to be likely, as mass motions 
have been detected even in quiescent active regions (Acton 
et al., 1981). We therefore attempted to evaluate the non- 
thermal terms as well as the thermal energy throughout the 
gradual phase of the five prime flares. 

As the flare observations made with SMM were not spe- 
cifically designed to address the energetics problem, our 
choice of flares (see Appendix 5 A) represents a compromise. 
We identified a number of different “types” of flares that, 
although sharing the same physical processes, had somewhat 
different physical characteristics. 

The Skylab Solar Flares Workshop (Sturrock 1980a) laid 
the foundations for the methodology used in this section, and 


Our goal was to reexamine these problems with the data from 
SMM and other supporting instruments as well as to take 
advantage of recent theoretical advances. SMM is capable 
of measuring coronal temperatures more accurately and with 
a better cadence than has been possible before. The SMM 
data set is also unique in that the complete transit of an ac- 
tive region was observed, with soft X-ray and UV images 
being taken every few minutes. We are therefore able to es- 
tablish the pre-flare conditions of the region and see whether 
anything has changed as a result of the flare. 

In the next subsection we describe the assumptions made 
in attempting to determine the required plasma parameters. 
The derived parameters for the five prime flares are pre- 
sented, and the role of numerical simulations is discussed. 
Finally, we consider the overall implications of our results 
and discuss how both theory and observations have evolved 
since the Skylab workshop. 

5.3.2 The Basic Physical Expressions 

The quantities needed for this study are defined in terms 
of the fundamental plasma parameters of electron tempera- 
ture (T e ), electron density (N e ), plasma volume (V), plasma 
velocity (v), and height (h). The four basic energies that we 
wish to obtain throughout the flares are the thermal energy 
(Eth) of the plasma, 

Et h = 3 N e k T e V ergs; (5.3.1) 

the kinetic energy (E k ) of the plasma, 

E k = 1/2 M p N e V v 2 ergs; (5.3.2) 

the potential energy (E p ) of the coronal plasma, 

Ep = M p N e V g 0 h ergs (5.3.3) 

(where M p is the proton mass and g 0 is the gravitational ac- 
celeration at the solar surface); and the stored ionization 
energy in the plasma (Ej), which is the sum of the ioniza- 
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tion energy of every ion stage of all abundant elements ex- 
isting up to the observed temperature. This latter component 
is dominated by the ionization energy of hydrogen and he- 
lium. Note that these energetics terms only apply in the above 
form for an isothermal plasma. However, the SMM data in- 
dicate that the coronal plasma during a flare is multi-thermal, 
and the integral forms of these equations must often be used 
(Section 5.3.2. 1). 

The energy flux terms and their corresponding timescales 
-are as important as the individual energy terms for under- 
standing the energy budget. There are three main terms: the 
conductive flux (F c ), which, in terms of the classical Spitzer- 
Harm theory of heat conduction, can be expressed as 


F c = — x T| /2 dT e /dh ergs cm -2 s -1 ; 

(5.3.4) 

the radiative flux (F r ), 


F r = — A~‘ j N| P r dV ergs cm -2 s _1 ; 

(5.3.5) 

and the enthalpy flux (F e ), 


F e = (5/2A) K T e d(N e V)/dt ergs cm 2 s -'. 

(5.3.6) 

The corresponding decay timescales (in s) are 


t c =2 X 10-'° N e L 2 T e 5/2 

(5.3.7) 

r r =10 7 T e /N e 

(5.3.8) 

T e — 1.2 N e (dN e /dt)-*. 

(5.3.9) 


It should be possible to derive all these quantities as a func- 
tion of time throughout the flare from a knowledge of the 
plasma parameters and the geometry of the flaring loop. 
However, these quantities are not available directly from ob- 
servations although the temperature and density can be in- 
ferred from soft X-ray line ratios. Generally, we derive the 
emission measure (N| V) from optically-thin line intensities. 
The SMM X-ray data can be used to derive the DEM. The 
temperature range available was 2 x 10 6 to 10 8 K, but no 
information was available in one of the most interesting tem- 
perature domains, namely 2 X 10 5 to 2 X 10 6 K. There is 
considerable uncertainty in the quantities needed to derive 
these energetics terms, especially the density. However, for 
the five prime flares, the values obtained for each of these 
energetics terms are discussed. 

5.3.2. 1 Multi-thermal Modeling 

Several quantities are discussed above which express the 
energy content of the various flare components with the as- 
sumption of an isothermal plasma. However, there is strong 
evidence from SMM data that the coronal plasma is multi- 
thermal during flares (Figure 5.3.1). For this reason, an at- 


tempt was made to model the temperature distribution of the 
flare plasma from spectroscopic data obtained with various 
SMM instruments. In this section we present the multither- 
mal forms of the expressions for the thermal energy (Eth) 
and radiative flux (F r ) and discuss the techniques which were 
used to derive the temperature distribution of the plasma. 
The temperature distribution of the plasma may be charac- 
terized by the differential emission measure (DEM), q(T), 
which is defined as 

q(T) dT = N|dV. (5.3.10) 

With a multithermal assumption, several of the quantities 
given in Section 5.3.2 are modified: the expression for the 
radiative flux becomes 

F r = J P r (T)q(T)dT ergs s ’ 1 , (5.3.11) 

where P r (T) is the radiative power loss function; the ther- 
mal energy, dE^, of a volume element dV with tempera- 
ture T and density N e is given by 

dEth = 3 N e kT d V ergs. (5.3.12) 

If the density within the volume, V, is assumed to be uni- 
form, Equation 5.3.12 may be integrated to yield 

E,h = 3kV l/z jTq(T)dT 

(Jq(T)dT) ,/2 efgS ' (5.3.13) 

Alternatively, if the pressure, p = nT, within the volume 
is assumed to be uniform, the thermal energy becomes, 

Eth = 3 k V" 2 [JT 2 q(T) dT]* /2 ergs. (5.3.14) 

The DEM was derived for the soft X-ray emitting ther- 
mal plasma by two different methods. The first method was 
an inversion technique initially developed by Withbroe (1976) 
and extended by Sylwester etal . (1980, 1984). This method 
has been applied to the five prime flares in this study by using 
data from four BCS lines and the six HXIS energy bands. 

The four lines observed by BCS consist of the helium- 
like resonance line OSq-Ti) in Ca XIX and Fe XXV and 
the lithium-like satellite line k in Ca XVm ( 2 P 1/2 - 2 D 3/2 ) 
and j in Fe XXTV ( 2 P3/2~ 2 D5/ 2 ). These satellite lines are 
dielectronically populated and so their intensity ratios rela- 
tive to the resonance lines are temperature sensitive and inde- 
pendent of ionization balance. The effects of blending with 
nearby lines were removed by fitting a fully synthesized spec- 
trum to the data. 

The HXIS coarse field-of-view data were analyzed by 
summing over the whole flare area rather than by studying 
individual pixels or groups of pixels, because it covers an 


5-21 


area comparable with that of the BCS collimator (6 arc min 
FWHM). This approach probably does not affect estimates 
of the overall energy budget of the flare but it could pos- 
sibly mask secondary energy releases at other sites away from 
the original flare location. 

Any particular solution of q(T) is probably not unique. 
However, various tests were performed which showed that 
different but acceptable, solutions for q(T) generally con- 
served the total emission measure, jq(T)dT= jN e dV. A fur- 
ther limitation arises from the fact that BCS and HXIS are 
insensitive to emission from plasmas with temperatures be- 
low about 10 7 K, although a considerable amount of thermal 
energy is expected to be contained in these plasmas. Both 
of these problems could be lessened by including softer X-ray 
or UV data, which would extend the results to lower tem- 
peratures and thus better constrain the range of acceptable 
solutions for the DEM. The lowest energy FCS channels are 
sensitive down to 2 x 10 6 K. However, since the cadence 
of the FCS observations was typically about 300 s per im- 
age, the FCS data were included for only one flare; 1980 
August 31 double flare. 

A second method for determining the DEM which com- 
bines data from BCS, HXIS, and HXRBS, has been de- 
veloped by Bely-Dubau et al. (1984). In contrast to the 
approach just described, this method assumes a parametric 
description for the DEM. The parameters are adjusted until 
an acceptable agreement with the observations is achieved. 
So far, this method has only been applied to the 1980 June 
29 flare. This technique is not subject to the problems which 
arise with an inversion technique when the emission func- 
tions cover a large proportion of the temperature range of 
the model. Eight BCS lines (the Ca XIX w, q, and k lines, 
Fe XXV w, q, and j lines, and Fe XXV r and s lines), the 
six HXIS energy bands, and the eight lowest energy HXRBS 
channels are included in the model. The model includes ther- 
mal and non-thermal components, and thus an important fac- 
tor for energetics studies is the lower-energy cut-off for the 
nonthermal spectrum. Bely-Dubau et al. (1984) found that 
any cutoff value between 15 and 25 keV gave an acceptable 
fit to the data, but resulted in a considerable variation in the 
estimates for the energy budget of the flare. 

Bommann (1985a, b) has developed a new method for es- 
timating the temperature and emission measure during the 
decay of the gradual phase from the shape of the light curves 
of soft X-ray lines. The rates at which the observed line fluxes 
decayed were not constant. For all but the highest tempera- 
ture lines observed, the rate changed abruptly, causing the 
fluxes to fall at a more rapid rate later in the flare decay. 
These changes occurred at earlier times for lines formed at 
higher temperatures. Bornmann proposed that this behavior 
was due to the decreasing temperature of the flare plasma 
tracking the rise and subsequent fall of each line emissivity 
function. This explanation was used to empirically model the 
observed light curves and to estimate the temperature and 


the change in emission measure of the plasma as a function 
of time during the decay phase of the 1980 November 5 flare. 
This method provides a simple and independent measure of 
the temperature and cooling rate during the decay of the flare. 

5.3.3 The Prime Flares 

A more detailed description of the prime flares and a list 
of publications discussing their relevant aspects are given 
in Appendix 5A. In keeping with the impulsive phase sec- 
tion, we have emphasized the 1980 April 8 flare, but a brief 
discussion of the other prime events is given below. 

5.3.3. 1 The 1980 April 8 Flare 

The M4 flare of 1980 April 8 originated in a small, iso- 
lated, bipolar region that was part of a delta configuration 
in NOAA active region 2372. The flare had an optical im- 
portance of IB and showed multiple impulsive spikes in the 
HXRBS light curve, the brightest occurring at 03:05 UT. 
For details see Appendix 5 A. 

Both BCS and FCS showed that until 10 min before the 
flare, the level of coronal emission from the whole region 
was declining (25 % in Ca XIX and about a factor of 2-3 in 
O VIE). The coronal temperature also appeared to be slowly 
dropping, changing from about 2.5 X 10 6 to less than 
2 X 10 6 K during the interval from 00:50 UT to 02:59 UT. 
At about 02:50 UT, the BCS Ca XIX lines started to increase 
in intensity. At 02:54 UT, the density measured with UVSP 
at the eastern footpoint began to rise exponentially, reach- 
ing a value of about 3 x 10 12 cm -3 at 03:02 UT, when the 
measurement sequence was terminated. Both intensity and 
density were still rapidly increasing at this time, so that con- 
siderably larger values could have been reached at the time 
of the peak emission. Machado (1980) has identified the 
eastern footpoint with the isolated bipolar region discussed 
above. The electron density, as derived from the P78-1 O 
VII line ratios, reached its maximum of 10 12 cm -3 at 03:05 
UT, before any of the soft X-ray intensities had peaked. Since 
the characteristic temperatures of the two determinations 
differ by an order of magnitude, the fact that the derived 
densities are comparable suggests that the two regions were 
not in pressure equilibrium with each other and implies the 
existence of temperature gradients across the magnetic field. 
This result must be treated with caution, however, since the 
two determinations were not made simultaneously. 

Subsequent development of the flare first involved the 
leader spot, via loops connecting it to the intermediate bi- 
polar region, after which it spread to the trailer spot. Later 
in the gradual phase, the region of soft X-ray emission ex- 
panded southward. The HXIS images show that at least three 
discrete sources were involved during the rise phase, sug- 
gesting the existence of at least two flaring loops. A descrip- 
tion of the magnetic field structure can be found in Chapter 
1 of these proceedings. 
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The observations were analyzed along the lines discussed 
in Section 5.2.3. The HXRBS data can be represented by 
either a thermal or a power law spectral form. For analysis 
of the microwave data from Toyokawa, temperatures deter- 
mined from fits to the thermal form of the HXRBS spectra 
were used to characterize the electron source of the hard 
X-ray emission for selected intervals throughout each burst. 
Observations obtained at 3.75 GHz were in the optically thick 
portion of the microwave spectrum throughout the impul- 
sive and gradual phases of the flare and therefore enabled 
the computation of an effective area. The 3/2 power of an 
effective area was taken as an estimate of the size of the emit- 
ting volume. The maximum turbulent velocity derived from 
the BCS spectra was 120 km s _1 , assuming a single 
temperature determined from spectral analysis and neglect- 
ing the source size (Antonucci etal ., 1982). The upflow was 
200-300 km s* 1 . Non-thermal line widths for the Fe XXV 
lines from P78-1 were derived by U. Feldman and S. Graham 
in work done specifically for this workshop. They assumed 
that the line profile was given by the convolution of three 
gaussians: one for the thermal width of the line, one for the 
turbulent velocity, and one for the angular size of the source. 
They assumed the angular size to be 30 arc sec and used the 
temperature value derived from P78-1 measurements of the 
j/w line ratio. Their results range from 100 to 300 km s 1 , 
bracketing the result of Antonucci et al . (1982). 

We used the DEM distribution to represent the ensem- 
ble of soft X-ray observations of the gradual phase of the 
flare. We chose this approach because the DEM is propor- 
tional to both the observed intensities and the total radiated 
power, even in the presence of inhomogenities in the radi- 
ating atmosphere. Four computations (Figure 5.3.1) of the 
DEM were made, corresponding to the four FCS rasters. 
BCS and HXIS data were selected for the times when the 
flaring plasma was being scanned with FCS. This technique 
is particularly powerful in placing constraints on the DEM 
distribution. The use of several isolated spectral lines pro- 
vides intensity measurements for which the contribution fane 
tion is well known. The inclusion of the four broad-band 
measurements from HXIS ensures that no important radia- 
tive energy sources have been omitted and provides a con- 
straint on the total emission. Each DEM distribution shows 
two peaks, one near 2 x 10 7 K and the other near 4 X 10 6 K. 
The high-temperature peak is seen to move toward lower 
temperatures as the flare cools, as one would expect. 

The temperature curves (see Figure 5.3.2) should be in- 
terpreted as indicators of the slope of the electron energy 
distribution in the energy range appropriate to each instru- 
ment; they do not necessarily represent physically meaningful 
quantities. With the exception of the temperature derived 
from HXRBS, the other temperatures agree to within a fac- 
tor of 2. Assuming these estimates are a measure of the 
plasma temperature, they are not misleading, provided that 
one does not attempt to extrapolate outside the relevant elec- 


tron energy range. The higher values of the HXRBS tem- 
perature parameter show that the electron spectrum producing 
the hard X-rays is relatively flatter than that producing the 
emission observed with the other instruments. 

Density measurements for this flare are available from 
several observational sources (Figure 5.3.3). At the time of 
the maximum at 03:05 UT, the O VII density is probably 
dominated by material in the footpoints and does not reflect 
the average density in the loop system, since the evapora- 
tion phase has just begun. 

Three volume estimates for the April 8 flare have been 
made: (1) the observed size in the HXIS fine field of view; 
(2) the apparent area derived from hard X-ray and micro- 
wave fluxes, and (3) using observed X-ray line intensities. 
The first two methods involve simple assumptions about the 
geometry of the emitting region to estimate the volume from 
the apparent area, whereas the third uses the derived elec- 
tron density to compute the effective volume from the ob- 
served emission measure. The two estimates based on 
apparent areas agree fairly well during the impulsive phase, 
whereas the volume calculated from derived density and 
emission measure estimates is lower by a factor of 100 to 
1000 (Figure 5.3.4). Although some of this discrepancy may 
be due to errors in the atomic data or in the underlying as- 
sumptions, the differences are important and have been cited 
as primary evidence for the existence of a small filling fac- 
tor. The volumes derived from all these methods converge 
during the gradual phase. This probably indicates evolution 
of the geometry of the plasma during the gradual phase, with 
the filling factor tending toward unity. The evaluation of the 
filling factor is of obvious importance to the question of ener- 
getics, since the thermal energy content of the plasma is 
directly proportional to the total volume. 


Figure 5.3.5 shows emission measures as functions of 
time determined from the BCS and HXRBS spectra. The 
P78-1 results of Doschek et al. (1981) are comparable. The 
four circles indicate the emission measures determined from 
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respond in time to the four FCS rasters. 

Values of P r and E^ (Figure 5.3.6) were determined for 
the high-temperature component (10 7 to 5 X 10 7 K) in a simi- 
lar manner as for the 1980 June 29 flare (see below). The 
volume of 10 27 cm 3 was estimated from HXIS fine-field-of- 
view images (Table 5.2. 1). The peak value for E^ was found 
to be about 1.4 X 10 30 ergs, which compares well with 
9.4 x 10 30 ergs for the total electron energy (Table 5.2.1). 
Note from Figure 5.3.1 that the low temperature component 
of the DEM is at least a factor of 10 larger than the high- 
temperature component which would mean that the thermal 
energy of the two components are approximately equivalent. 

Figure 5.3.7 is the version of Figure 5.1.1 for this par- 
ticular flare and shows the complete energy budget derived 
from the above physical quantities plotted. The kinetic energy 
in the upflows and in the turbulence was calculated using 
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COMPOSITE DIFFERENTIAL EMISSION MEASURE 
1980 APRIL 8 



ELECTRON TEMPERATURE (10 6 K) 


Figure 5.3.1 The DEM models calculated for four times during the 1980 April 8 flare. Note that, even during 
the rising portion of the gradual phase when the flare is at its hottest, there is a significant low temperature 
component. As the flare evolves the high temperature component of the DEM decays in temperature eventually 
merging with the low temperature component. 


the P78-1 velocities and the BCS emission measures. The 
energy in the electrons determined from the thermal fit to 
the HXRBS spectra was used in Figure 5.3.7 rather than the 
thick-target non-thermal energy since the instantaneous 
energy in the electrons for that case is not significant. The 
gravitational potential energy and ionization energy are not 
plotted here, as they are too small to show on this scale. The 
peak gravitational energy is 3 X 10 28 ergs. 

Note that at 03:06 UT the total energy remains essen- 
tially constant. The slow drop in the total energy after this 
time (less than 30%) may be due to conductive losses which 
have not been included in the plot. 


5.3.3.2 The 1980 May 21 Flare 

This flare was chosen because it has exceptionally good 
data for the impulsive phase. However, the FCS and UVSP 
coverage during the gradual phase was very poor. They were 
in unsuitable modes during this event and so can contribute 
very little to the discussion of the energetics of the gradual 
phase, which is well described by Duijveman et aL (1983). 

5.3.3.3 The 1980 June 29 Flare 

This event occurred on the limb, and so the determina- 
tion of the volume of the flare is different than for a disc 


5-24 







TIME (UT> 

Figure 5.3.2 Evolution of the temperature parameter for 1 980 
April 8 derived from a number of different sources. The dia- 
monds are derived from the high temperature component of 
the DEM calculations shown in Figure 5.3.1. The error bars 
represent the full width at half maximum of the distribution. 



Figure 5.3.3 Evolution of the density parameter from 
UVSP Si IV/O IV line ratios and O VII measurements 
from P78-1 , compared to the density derived from the 
hard X-ray — microwave analysis. 



Figure 5.3.4 The history of the volume parameter 
from the 1980 April 8 flare from several different 
sources. 
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TIME (UT) 

Figure 5.3.5 The evolution of the emission measure 
in the 1980 April 8 flare from the BCS, HXIS, and 
HXRBS data (thermal assumption). Note that the soft 
X-ray emission measure is over 3 orders of magnitude 
larger than that for the hard X-ray source. The BCS 
and HXIS value does not include information on the 
low temperature component (< 10 7 K). The filled cir- 
cles represent the emission measure calculated for 
temperature >2 x 10® K by including the FCS data 
with the BCS and HXIS observations. They show that 
the component between 2 x 1 0® K and 1 0 7 K is com- 
parable to the component in the 10 to 20 x 10 6 K 
range, at least later in the flare. 



Figure 5.3.6 The evolution of (a) the radiated power 
(ergs s _1 ) and (b) the thermal energy (ergs) for the 
1 980 April 8 flare from combined BCS and HXIS data. 
The error bars represent ± 1o uncertainties. These 
quantities are calculated for the thermal plasma above 
10 7 K only. 


Table 5.3.1 Energetics of the 1980 June 29 Flare 


Energetics term 

BCS + HXIS 

Ca XIX 
Fe XXV 

Fe XXI 
0 V 

HXIS 

HXRBS 

Emission measure 
(cm" 3 ) 

2 x 10 49 

3.5 x 10 49 

5.6 x 10 49 

7.9 X 10 43 
0.32 x 10 43 

3.3 X 10 45 

9 x 10 45 

Radiative power 

8 x 10 26 

2 x 10 27 


2.3 x 10 28 



(ergs s" 1 ) 

Integrated thermal 2.8 x 10 26 

energy (ergs) 

Radiated energy 3.7 x 10 29 

(ergs) 


5-26 




10 ' 


i r 

-_m RADIATION 
lo THERMAL PLASMA 
- m ELECTRONS 
UPFLOWS 


I 

EH TURBULENCE 
/77 ACTIVE REGION 
THERMAL PLASMA 



03:05 


03:10 

TIME (UT) 


03:15 03:20 


Figure 5.3.7 The distribution of the flare energy 
among various reservoirs and active phenomena as 
a function of time for the 1 980 April 8 flare (c.f . Figure 
5.1.1). The width of the histogram bins represents the 
time resolution for each of the parameters measured. 
Note that late in the flare (after 03:09 UT), the sum 
of the components remains effectively constant in 
spite of the fact that several of the components are 
varying. 


flare. Although we have the height structure and hence the 
area of the flare, the line-of-sight depth is still unknown. Fur- 
ther, it seems that part of the flare was probably behind the 
limb. The FCS images were used to determine the area of 
the flare, which gave a value of 2 x 10 19 cm 2 . The micro- 
wave data imply an area of 3.3 x 10 19 cm 2 , although inter- 
ferometric observations gave an area of an order of 
magnitude less. 

Normally, one does not expect to find significant line shift 
in the BCS data for a limb event. However, this flare gave 
a velocity of 370 km s -1 at the flare onset. Turbulent veloc- 
ities of about 200 kms" 1 were estimated from the BCS Ca 
XIX line profiles. Vertical motions were measured by several 
instruments: 200 km s' 1 in the lower corona, 500 km s _1 
above 1.2 solar radii, and 1100 km s _l from the type II 
burst. 


A multi-thermal analysis of the high-temperature com- 
ponent of the flare plasma was made for this flare. For the 
determination of the DEM, the intensities of four lines ob- 
served by BCS and the integrated HXIS coarse-field-of-view 
fluxes were used. The DEM was then convolved with the 
radiative power loss function, P r (T), of Summers and 
McWhirter (1979) to give the total radiative loss rate results 
shown in Figure 5.3.8(a). Since P r (T) is approximately con- 
stant with temperature near 2 x 10 7 K, the behavior of the 
radiative flux with time is similar to the behavior of the total 
emission measure with time. The value of E^ is shown in 
Figure 5.3.8(b), as calculated assuming either uniform pres- 
sure or uniform density throughout the volume of plasma 
with a given temperature. The volume was estimated to be 
7 X 10 26 cm 3 by assuming a density of 2 X 10 11 cm 3 , 
which is typical of the estimated densities for the other prin- 
cipal flares we considered (see Table 5.2.4). The peak value 
of 7.8 x 10 29 ergs for E^ calculated assuming uniform den- 
sity is comparable to the total energy inferred (1.3 X10 30 
ergs) for the non-thermal electrons above 25 keV. The 
turbulent energy derived from the BCS line widths was dis- 
cussed in Section 5. 2. 3.4. The peak turbulent energy ob- 
served during the impulsive phase was approximately 
3 X 10 28 ergs. 

There was a mass ejection during the 1980 June 29 event 
(Wu et al ., 1983). The maximum temperature (2 X 10 6 K) 
and density (4 X 10 u cm~ 3 ) derived from the XRP data, 
together with an assumed velocity, were used to simulate 
the mass ejection using an MHD model (Wu et al ., 1983). 
The simulation produced spatially wide, large-amplitude, 
temporally steepened MHD waves but no shocks. This result 
seems to be supported by the fact that a type II radio burst 
was observed late in the event for only a few minutes. The 
density enhancements moved away from the Sun at the same 
velocity (1000 km s _1 ) as determined with the Mauna Loa 
K-coronameter. However, the observation of the coronal 
transient included a rarefaction that does not appear in the 
simulation. The model enables us to derive the energy as- 
sociated with the transient; the quantities are: peak temper- 
ature, 16.5 X 10 6 K; peak density, 3.0 X 10 10 cm' 3 ; total 
volume, 1.0 x 10 27 cm 3 ; and total energy, 4.0 X 10 33 ergs. 
Comparison with the observed pB (polarization X bright- 
ness) signal, derived from C/P data to be 7 X 10 8 , gave 
a result that was a factor of 3 too high. This implies that 
the magnitude of the calculated density enhancement was high 
by at least 50% in the closed-field case and a factor of 3 
higher in the open-field case. The peak values of the ener- 
getics terms of this flare are listed in Table 5.3. 1 , in which 
only part of the radiative energy could be estimated based 
on the available data. The computed total energy included 
all the energy sources (i.e. magnetic, potential, kinetic and 
thermal). 

Using the alternative technique of Bely-Dubau et al . , 
(1984) described above (Section 5.3.2. 1), we find the energy 
content as a function of time given in Table 5.3.2. 
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Figure 5.3.8 The radiated power (a), thermal energy 
(b), and kinetic energy (c) of the soft X-ray emitting 
plasma as functions of time for the 1980 June 29 
flare. These quantities were derived from the BCS and 
HXIS data. The error bars represent ± 1o uncertain- 
ties. The error bars on the kinetic energy estimates 
were calculated from the uncertainties of the line 
widths and electron temperatures but neglect the un- 
certainties of the estimated volume and emission 
measures. 


S.3.3.4 The 1980 August 31 Flare 

The properties of the impulsive double flare on 1980 
August 31 were reported by Strong et al. (1984). By com- 
bining observations of the two events from several of the 
SMM instruments (XRP, HXIS, HXRBS, and UVSP) as well 
as radio and optical observatories, they concluded that the 
second energy release probably occurred in the corona 
plasma that was hotter and denser than that of the first flare. 
They inferred that the differences between the two flares 
resulted from changes in plasma properties brought about 
by the first energy release. Comparisons of the radio and 
X-ray data show that the primary energy release site was con- 
tiguous (for energetic electrons) with separate, closed mag- 
netic structures. The loops which emitted the radio bursts 
were at least 10 times longer than those emitting the X-rays. 
There was no evidence for emerging magnetic flux in this 
region around the time of the flare; in fact, the field seemed 
to be decaying rapidly as the two leader spots approached 
each other. The primary energy release mechanism involved 
therefore seemed to be magnetic reconnection caused by the 
“collision” of these structures due to photospheric motions. 

The density was determined from the emission measure 
and source size. The pre-flare temperature of the active 
region was found to be approximately 2.5 X 10 6 K with an 
average emission measure of 3 x 10 47 cm -3 . As there were 
no distinct structures visible, the volume was derived from 
the area of the FCS pixel and a typical soft X-ray scale height 
and was found to be 4 x 10 28 cm 3 . The mean electron den- 
sity was therefore 3 X 10 9 cm -3 for the pre-flare active 
region. Although there were no pre-flare magnetograms, data 
taken just after this event showed that there was ample free 
magnetic energy available even after the flares, to supply 
all the energy observed. Over the previous days, observa- 
tions of this active region had indicated that the two leading 
spots (see Section 5A.4) were moving towards one another. 
Although the velocities were relatively low (of the order of 
hundreds of meters per second), the energy involved in these 
photospheric motions must have been large and supplied the 
stored magnetic energy in the form of sheared magnetic 
fields. 

The highest spatial resolution available on SMM (3 arc- 
sec from UVSP), indicated that the X-ray emitting region 
in the first flare was 10 26 cm 3 with an uncertainty of about 
a factor of 2. The same volume was assumed throughout both 
flares, as there was no direct measure of the volume for the 
second flare although the longer timescales of the second flare 
imply an increased volume. 

The multi-thermal analysis covered the extended temper- 
ature range from 2 X 10 6 K to 5 X 10 7 K by using six FCS 
channels, four BCS lines, and the six HXIS fme-field-of-view 
channels. Figure 5.3.9 shows P r , E th and E^ as functions 
of time. The thermal analysis was not possible until the time 
of the peak of the first flare in soft X-rays as the count rates 
were low. In calculating E th , a volume of 7.6 X 10 25 cm 3 
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Table 5.3.2 Energy as a Function of Time for the 1980 June 29 Flare 


Time (UT) 

18:20:40 

18:23:28 

18:23:54 

18:25:00 

18:26:00 

18:27:18 

Thermal 







Energy 
(10 30 ergs) 

1.38 

3.4 

5.5 

7.6 

7.6 

8.3 

Energy 
Deposit rate 
(10 27 ergs s _1 ) 

1700 

2000 

26 

30 

4 

2.5 

Cutoff 

energy (keV) 

10 

12 

26 

26 

26 

26 



Figure 5.3.9 Same as Figure 5.3.8 for the two flares on 1980 August 31. Compare the peak kinetic energy 
from the second flare with that of the second flare on 29 June. There is nearly an order of magnitude difference 
showing that the role of mass motions in this second event are very much inhibited compared to both that 
of the first flare and the flare on 29 June. 
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estimated from the HXIS fine-field-of-view images was used 
for both flares (see Table 5.2.2) even though the volume of 
the first flare was probably three times larger. The peak value 
of Eth for the second flare was approximately 1.6 X 10 29 
ergs and compares with a total electron energy of 5 X 10 29 
ergs. 

The energy associated with the turbulent broadening was 
discussed in Section 5. 2. 3.4. The peak turbulent energy dur- 
ing the second flare was about 6 x 10 27 ergs. It is interest- 
ing to note that the ratio of the peak value of to the peak 
turbulent energy for this flare is approximately 27, which 
compares with 26 for the same ratio for the 1980 June 29 
flare. This ratio is close to 10 for the first flare on 1980 
August 3 1 . 

A comprehensive picture of the mass motions was derived 
from the BCS X-ray line profiles and from the source mo- 
tions observed by HXIS. The uncertainties in the velocity 
measurements are minimal because the flare was compact 
and therefore not complicated by multiple flare sites. Using 
the peak integrated DEM and the above volume of 7.6 X 
10 25 cm 3 , we estimated the mass of the “evaporated” plasma 
and hence the kinetic energy. The potential energy was de- 
termined at the peak of the soft X-ray emission from the 
source size and assuming a semicircular geometry for the 
loops. Again, by using the density and temperature, we es- 
timated the ionization energy. The results of these calcula- 
tions are listed in Table 5.3.3 for both flares. 

This event illustrates multiple energy release, seemingly 
a common aspect of flares and a possible reason for so much 
diversity in flare characteristics. It establishes, for some cases 
at least, the corona as the locale for the primary energy 
release involving closed magnetic structures. It also shows 
that structures of vastly different scale sizes can be involved 
in the same impulsive release. These are the sorts of morpho- 


logical boundary conditions which are a prerequisite for con- 
tinued progress in flare theory. 

5.3.3.5 The 1980 November 5 Flare 

This was one of the most comprehensively observed 
flares, with good groundbased data as well as good configu- 
rations for the SMM instruments. It was another example 
of a double flare but with the components separated by more 
time than for the 1980 August 31 events. The geometry of 
the flare has been derived by Martens et al . (1984), who 
found the two principal loops observed during the flare to 
be 1 .6 x 10 9 cm and 7 x 10 9 cm in length, with a diameter 
of 1 . 1 x 10 9 cm. The volume of the larger of the loops was 
therefore about 4 x 10 28 cm 3 if it is assumed to be semi- 
circular. 

Using the higher density values from the He-like ion ratios 
(discussed below), we derived the energetics terms shown 
in Table 5.3.4. Figure 5.3.10 shows P r and E^ for the two 
flares on November 5. Estimates of the temperature and 
emission measure were made from the following parameters: 
the shape of the FCS soft X-ray light curves (Bornmann 
1985a,b), ratios of GOES count rates, DEM fits to the FCS 
data, satellite-to-resonance line ratios observed with BCS, 
ratios of HXIS count rates, and thermal fits to HXRBS data. 
The temperature rose to about 4 x 10 8 K in both flares and 
fell rapidly during the initial decay, relaxing to a nearly ex- 
ponential decay late in the decay phase. The emission meas- 
ure derived from HXRBS reached a maximum of 8 X 10 45 
cm -3 and 10 47 cm -3 for the first and second flare, respec- 
tively. The soft X-ray instruments showed emission meas- 
ures which began and ended near 10 48 cm -3 and reached a 
maximum value of 3 to 6 x 10 49 cm -3 . 

The electron density estimates during the second flare 
differ by roughly two orders of magnitude. Duijveman et 


Table 5.3.3 Non-Thermal Terms of the Energy Budget for the 1980 August 31 Flares 


Energy term 

Kinetic energy 
Turbulence 
Upflows 
Translational 
Total 


Energy in first flare (ergs) 

4 x 10 29 

5 x 10 28 
7 x 10 28 
5 x 10 29 
1 x 10 27 

3 x 10 26 


Energy in second flare (ergs) 

< 3 x 10 27 

< 2 x 10 27 
1 x 10 28 
1 x 10 28 

1.5 x 10 27 


Potential energy 

Ionization 
potential energy 
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3 x 10 26 


Table 5.3.4 Energetics Components in ergs 
for the 1980 November 5 Flare 


Term 


Peak Decay Total 


Radiative losses 3 x 10 28 

Thermal energy 7 x 10 29 

Mass motions (upflows) 4 x 10 29 

Mass motions (turbulence) 4 x 10 28 


6.2 X 10 26 
3 x 10 27 
2 x 10 28 


2.8 x 10 29 



TIME (UT) 

Figure 5.3. 1 0 The evolution of (a) the radiated power 
(ergs s _1 ) and (b) the thermal energy (ergs) for the 
1 980 November 5 flare from combined BCS and HXIS 
data. The error bars represent ± la uncertainties. 
These quantities are calculated for the thermal plasma 
above 10 7 K only. 

al. (1982, 1983a) used the emission-measure volume tech- 
nique, assuming a line-of-sight depth of the emitting region 
of 6000 km, to find electron densities in the two principal 
loops of 2.5 X 10 10 and 9 X 10 9 cm -3 . Density estimates 
made from FCS line ratios (Figure 5.3.11) indicate that the 
density fell from 1.5 x 10 12 cm -3 at the time of peak soft 
X-ray flux to 3.0 x 10 u cm -3 near the end of the flare 
decay (Wolfson et al., 1982; Bornmann 1985a). For times 
when both density estimates are available, the emission- 
measure volume estimates are a factor of 30-50 lower than 



Figure 5.3.11 The evolution of the density and the 
Ne IX flux during the flare of 1 980 November 5. The 
density is derived from the line ratios of He-like Ne IX 
(figure courtesy of C.J. Wolfson). 


the line ratio estimates. This discrepancy could be the result 
of emission from separate spatial regions, but this seems un- 
likely since both detectors sampled overlapping regions in 
space and energy. Alternatively, the line ratio estimate 
could be correct if the filling factor was 10~ 3 to 10 -4 . 

However, Wolfson et al. (1983) found that the emitting 
volume observed with FCS decreased from 3 x 10 25 to 
3 X 10 24 cm 3 during the decay of the flare, assuming the 
temperature remained constant. The volume of an FCS pixel 
is about 10 27 cm 3 (A 3/2 ), which gives a filling factor of be- 
tween 3 X 10~ 3 and 3 x 10 2 . Thus, two separate esti- 
mates of the filling factor for this flare both indicate the need 
for filling factors of the order of 10~ 3 . 

DEM analysis following the method of Withbroe et al. 
(1975) and Sylwester et al. (1980) using the FCS data indi- 
cates a high-temperature component which diminished in 
magnitude and temperature during the decay phase. A sec- 
ond, low-temperature component remained nearly constant 
with time. This DEM behavior has been observed in other 
flares examined by our group (cf. Figure 5.3.1). 

Bornmann ( 1985a, b) used the shape of the FCS light 
curves to determine the temperature and emission measure 
during the decay of this flare (Figure 5.3.12). She assumed 
that the plasma was isothermal and that discrepencies be- 
tween observed and calculated fluxes (which were less than 
a factor of 2) were due to incorrect assumptions for detector 
efficiencies and/or atomic abundances. Although the DEM 
and light curve analyses both used the same data, they lead 
to different conclusions because of the manner in which the 
two techniques respond to possible problems in absolute flux 
calibrations; the DEM method assumes that the values given 
are correct, whereas the light curve method assumes that the 
plasma is isothermal and adjusts these values accordingly. 
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Figure 5.3.12 FCS fluxes for the 1980 November 5 flare illustrating the sharp changes in slope for each of 
the ions. Note how the ions formed at the higher temperatures have breaks earlier in time. 






The total emission measures in both cases agree to a factor 
of 2. The temperature derived from the light curve lies be- 
tween the temperatures of the two peaks in the DEM. Fur- 
ther study of these discrepancies is in progress. 

The volume derived from the soft X-ray emission meas- 
ure and electron density (determined from line ratios) was 
found to increase from 7.2 x 10 24 to3.4 x 10 25 cm 3 during 
the flare decay (Bommann 1985a). The spatial extent of the 
flare, as recorded in the lowest-energy HXIS channel (3.5 
to 5.5 keV) was not changing during this time, at least on 
the scale of the 8-arc-sec spatial resoluton. This indicates 
that the filling factor increased with time during the flare 
decay. The work done, if the increase in volume is due to 
expansion of the emitting material rather than to the addi- 
tion of emitting material in neighboring regions, was W = JP 
dV > jng k T dV, or greater than 3.2 x 10 28 ergs (Bom- 
mann 1985a). 

Bornmann (1985a) concluded that the mass of the soft 
X-ray emitting plasma remained constant during the decay 
phase, although the large errors on this measurement do not 
eliminate the possibility that the mass varied by a factor of 
two, reaching a maximum late in the decay phase. 

The total radiative loss rate (Nagai 1980, Rosner et al . , 
1978) for the equivalent isothermal plasma derived from the 
temperature and emission measure was found to fall from 
6.3 X 10 26 to 1.9 x 10 25 ergs s" 1 during the decay of the 
second flare (Bornmann 1985a). The total radiative loss dur- 
ing this period was 2.8 x 10 29 ergs. 

The total thermal energy (E^ = k T V) of the plasma 
seen in the FCS pixel during the decay of the second flare 
was 4.2 X 10 28 ergs at the time of peak soft X-ray flux, 
decreasing to 2.2 X 10 28 ergs (Bommann 1985a). Because 
FCS sampled only a small portion of the flare, the total ther- 
mal energy of the entire flare plasma was probably an order 
of magnitude larger. The thermal energy in the decay phase 
fell at a rate of 2.7 x 10 25 erg s _1 (Bommann 1985b), 
which is slower than the rate at which energy was lost by 
radiation. Therefore, additional heating during the decay 
phase is required. 

The energy in mass motions fell off rapidly, at the rate 
of 1.8 X 10 28 ergs s~ l . The conversion of the energy in 
mass motions into the observed radiative and thermal ener- 
gies could be the source of additional heating, but the rate 
of decay of the mass motions was too rapid and would re- 
quire an intermediate storage mechanism before appearing 
as thermal or radiant energy. The turbulent line broadening 
began at the same time that the blue shifts appeared and con- 
tinued for several minutes. The BCS line widths indicate that 
the turbulent velocities peaked at 100 to 200 kms" 1 , simul- 
taneously with the hard X-ray peak, and fell essentially to 
zero during the decay phase. The peak turbulent energy, as- 
suming an electron density of 10 12 cm~ 3 , was 5.6 X 10 27 
ergs and fell at the slower rate of 5 x 10 25 ergs s' 1 . If this 
rate continued past the time at which turbulent broadening 


was observable, it could account for a significant fraction 
of the energy lost by radiation. 

5.3.4 Flare Modeling 

The above studies do not address the problem of the ener- 
gy source function. Numerical simulations require an energy 
input function to reproduce the observed fluxes in the coronal 
lines. By adjusting the input energy profile to agree with the 
observed evolution of the X-ray line intensities, it should be 
possible to derive the heating function. The recent rise in 
interest in numerical simulations of solar flares is dealt with 
in more detail in Chapter 7. However, that section empha- 
sizes the numerical techniques and differences between the 
codes. In contrast, the work done for our group by R. 
Pallavicini, G. Peres and S. Serio used a hydrocode to simu- 
late two real flares. 

The 1980 May 7 (15:00 UT) and November 12 (17:00 
UT) flares were chosen for this study. In both cases the ini- 
tial conditions for the models were determined from obser- 
vations. The one-dimensional, time-dependent fluid equations 
were solved (see Peres et al . , 1982; Pallavicini et al . , 1983) 
subject to Gaussian energy input of the form 

H = H 0 f t exp[-(s-so) 2 /(2 a 2 )]. (5.3.15) 

In the first stage of the analysis, sq was chosen so that the 
pulse was centered on the loop summit. In other work 
(Pallavicini et al., 1983; Peres et al., 1984), the fitting has 
been done for the pulse centered at the base of the loop and 
for energy deposition with electron beams. However, Pal- 
lavicini et al. (1983) showed that the exact location of the 
energy input did not make a significant difference to the 
global form of the computed light curves. The duration and 
amplitude of the heating during the rise phase of the flare 
are obtained by a best fit to the temperature and Fe XXV 
light curve. This approach is based on the results of Pal- 
lavicini et al., who showed that the Fe XXV light curves 
start to decay as soon as the heating level is decreased. 

The May 7 flare was a compact, short-lived event. FCS 
observations show that it was confined within a single 15 
arcsec pixel. An extensive analysis of this event, based on 
simultaneous SMM and Ha observations, has been given by 
Acton et al. (1982). The flare consisted of a single compact 
loop with a length (L) of 10 9 cm and cross-sectional area 
(A) of 9 X 10 16 cm 2 . The pre-flare temperature at the same 
location was determined to be T e = (2.5 ±0.5) x 10 6 K 
from the ratio of the intensites of the FCS Ne IX and Mg 
XI lines. The initial pressure and density have been inferred 
from the loop temperature and length by using the scaling 
law of Rosner et al. (1978). The derived initial pressure (Pq) 
was about 10 dynes cm -2 , implying a pre-flare coronal den- 
sity of 10 10 cm -3 . 

The parameters and the time evolution function used to 
fit the observations were Hq = 175 ergs cm' 3 s” 1 ; s 0 = 


C- v 
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5.7 X 10 8 cm; 5 = 10 8 cm; f(r) = t/150 for 0 < t < 150 s 
and f(t) = exp[(150 — t )/r] for t > 150 s. 

In the initial model the decay time (t) was assumed to 
be zero (i.e., an instantaneous switch-off of the heating func- 
tion). This allows a good fit to the BCS Ca XIX and Fe XXV 
light curves during the rise phase of this flare (Figure 5.3.13). 
The lower-temperature lines from FCS cannot be used to fur- 
ther constrain the fit, as the temporal resolution available 
during this time was only 155 s. These parameters do not 
fit either of the BCS lines during the decay phase. The 
predicted fluxes decay faster than the observed ones. The 
decay time was increased to accommodate this problem; a 
value of about 60 s would give very likely the desired result. 
The large increase in density and temperature in the flare 
generally causes a decrease in all the characteristic timescales 
for a given loop model (e.g. , the radiative cooling time, the 
conductive cooling time, and the sound travel time). This 
imposes a practical limitation on how far the evolution of 
the flare can be followed into the decay; in this work the 
model was limited to 20 s after the flare peak. 

The thermal energy derived from these models was, typi- 
cally, an order of magnitude larger than the gravitational and 
kinetic energies. From the modeling, the energy deposited 
during the rise phase of this flare was about 3 X 10 29 ergs. 
Allowing for an additional energy deposition during the de- 
cay phase with r = 60 s, the total energy input was 5 X 


10 29 ergs. The energy deposition would increase if the ener- 
gy were deposited at the loop footpOints. (For further de- 
tails see Peres et al (1985)). 

The November 12 flare was substantially different from 
the May 7 flare in that it had a larger spatial extent (cover- 
ing at least three 15-arc-sec pixels) and a somewhat lower 
density and pressure in the pre-flare state. For a more detailed 
description of this flare, see MacNiece et al. (1985). The 
soft X-ray light curves appear to be typical of most flares 
observed by FCS. Since the raster cadence was faster than 
for the May 7 flare, the simulations could be carried further. 

The flare occurred in a loop with L = 4 X 10 9 cm and 
A = 2.5 x 10 17 cm 2 . The pre-flare density and temperature 
were calculated as above to be 7 x 10 9 cm -3 and 3 X 
10 6 K, respectively, corresponding to a Pq of about 6 dynes 
cm -2 . The heat input was assumed to be at the loop sum- 
mit (hence so was 2 x 10 9 cm). The* input was assumed to 
be 10 ergs cm -3 s _1 at a steady rate for a period of 180 s 
with a value for s of 5 x 10 8 cm during the rise phase of 
the flare. Several values of r (0, 30, 60, and 140 s) were 
tried to fit the decay phase. The best fit to the observations 
was r = 60 s, which was mainly determined from fitting 
the Fe XXV light curve, owing to the strong temperature 
dependence of the emissivity of this line. 

Figure 5.3.14 shows a comparison of observations and 
numerical simulations for the FCS and BCS light curves. The 


7 MAY 1980 




SECONDS 

Figure 5.3.13 Predicted light curves (solid lines) from a loop model with an instantaneous switch off of the 
transient heating at t = 1 50s, compared to BCS data (crosses) for Ca XIX and Fe XXV (t = 0 corresponds 
to 1 4:56:36 UT). The disagreement between observed and predicted light curves in the decay phase indicates 
that additional heating is required after the peak in Fe XXV. 
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12 NOVEMBER 1980 




SECONDS 

12 NOVEMBER 1980 

Figure 5.3.14 Comparison of the predicted (solid lines) and observed (crosses) fluxes for the XRP soft X-ray 
lines in the 1 980 November 1 2 flare. The model has steady heating to t = 1 80s, followed by an exponential 
decay with an e-folding time of 60s (t = 0 corresponds to 1 7:00:00 UT). (a) Comparison for the six FCS chan- 
nels. (b) Comparison with the BCS lines. 
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predicted fluxes appear to be in agreement with observation 
for all lines except Ne IX. Not only is the general time evolu- 
tion well reproduced, but the absolute flux values also agree 
with observations within the estimated calibration accuracy. 
The discrepancy found for Ne IX is not readily explained 
except in terms of instrumental effects, such as blending with 
nearby Fe XIX lines and/or crystal fluorescence. Note that 
the choice of r = 60 s is made to fit only the Fe XXV light 
curve. 

From the above results, it is possible to evaluate the 
energy budget of the November 12 flare. Figure 5.3.15 
shows Eth as a function of time. Again, dominated all 
the other sources of energy by an order of magnitude at any 
time during the flare. The total energy deposited into the loop 
for r = 60 s was 7 X 10 29 ergs. 



Figure 5.3.15 The evolution of the thermal energy in 
the 1980 November 5 flare as predicted by the dy- 
namical loop model. 


The decay time we have adopted for the heating term has 
little direct relationship with the decay times of the light 
curves of the lines of lower ionization stages. In fact, the 
evolution of the light curves of the “cooler” ion is domi- 
nated by the evolution of emission measure and temperature 
(remember that during evolution, the loop’s temperature is 
hardly uniform). Instead, as soon as the heating is switched 
off, the Fe XXV radiative efficiency is a very steep func- 
tion of the temperature up to 70 x 10 6 K (see Figure 8 of 
Pallavicini et al., 1983). On the other hand, the emission 
measure of the cooler plasma (say at temperatures near the 
excitation temperature of O VIII) will keep increasing, es- 
sentially because the cooling of the hotter component will 
enrich the emission measure of the cooler one, if the maxi- 


mum temperature in the loop is less than 70 X 10 6 K. There- 
fore, a decay of the Fe XXV line is a good indication of a 
decay in the heating term, due to steepness of the Fe XXV 
radiative efficiency. 

Actually, the simulations for the November 12 flare, if 
the heating is switched off abruptly, Fe XXV decays much 
faster than observed. For this reason we need some smoother 
decay for the heating. An exponential decay with r = 60 s 
gives the best results for Fe. A posteriori, we also see that 
it gives agreement with the cooler lines, although for these 
last lines the value of r, as explained, is less significant. For 
the May 7 flare light curves, we can infer that r — 0 does 
not agree with the observations. The appropriateness of r 
= 60 s is only guessed by the shape of the light curves and 
by our experience of many other simulations. 

In summary, if one allows for heating in the decay phase 
of flares, good agreement with observations can be achieved. 
The success of this particular modeling effort lies in the fact 
that a reasonable fit to the form and, in most cases, the ampli- 
tude of the light curves of the remaining BCS and FCS lines 
was obtained, even though the parameters were selected to 
give a good fit only to the Fe XXV light curve. These results 
indicate that the bulk properties of the decay phase can be 
represented quite well by such a fluid model. 

5.3.5 Discussion 

So far in this section, we have outlined the methods of 
determining the physical parameters of our five prime flares 
and have discussed the energy budget at different times dur- 
ing those events. Having determined the temperature, den- 
sity, and length of the flare loops, we could follow Moore 
et al. (1980) and construct a diagnostic diagram based on 
the order of magnitude cooling times (Equations 5. 3. 7-5. 3. 9) 
to determine whether our flares require long-term heating. 
We follow a similar procedure, but first it is important to 
record the progress that has been made in our theoretical un- 
derstanding of flare cooling in recent years. 

An important physical property of the flare plasma is how 
the DEM scales with temperature. We find that for 5.5 < 
log(T e ) < 7 ( Dere and Cook, 1979; Widing and Spicer, 
1980; Acton et al ., 1983), 

q = A N e 2 T (dT/ds)-' * Tf, (5.3.16) 

where A is the emitting cross sectional area and f is typically 
3.5. If we can find a cooling mechanism that gives f = 3.5, 
we may be able to find out how flares cool. Antiochos and 
Sturrock (1976) considered static conductive cooling in which 
the downward heat flux is dissipated in the chromosphere, 
presumably by radiation. They also included a non-uniform 
magnetic field structure in which the field is stronger at the 
footpoint of the loop. Underwood et al. (1978) found that 
f = 0.5, and for the area change to have significant effect 
on the cooling time (a factor of about 3), the magnetic field 
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had to change by a factor of 20 from the loop top to bottom. 
A refinement of this model was made by Dowdy et al. 
(1983), who considered a variety of flux-tube geometries, 
and concluded that the shape and amount of divergence were 
important in determining the magnitude of the inhibition of 
the heat flux. A slightly different approach to conductive 
cooling was taken by Van Hoven (1979), who considered 
the static cooling in a twisted flux tube with constant cross- 
sectional area and showed that the more the flux tube was 
twisted, the more the conduction was inhibited. Clearly, 
twisting gives rise to longer field lines, which makes con- 
duction a less efficient mechanism. However, none of these 
variations in static conductive cooling makes any change in 
the form of q. Antiochos and Sturrock (1978) relaxed the 
static assumption and considered a model in which the chro- 
mosphere could respond to a downward heat flux by mov- 
ing (or evaporating) upward. This model differs in an 
important way from that in their previous paper (Antiochos 
and Sturrock 1976), since the thermal energy of the flare 
plasma is conserved, so that all the energy in the downward 
heat flux at the base of the model goes into driving an up- 
flow. They found that f = —0.5 but noted that no energy 
is lost in radiation either implicitly or explicitly, so that the 
effect of evaporation will be overestimated. They did not pro- 
vide a comparison between static and dynamic cooling. Since 
the conductive models did not appear to fit the slope of the 
DEM, Antiochos (1980) then tried radiative models. He 
found values of f between 1 and 2, depending on the form 
of the loss function chosen. Some of these models include 
a subsonic downflow, and as an extension to this, Antiochos 
and Sturrock (1982) included steady-state supersonic down- 
flows driven by a pressure gradient between the loop sum- 
mit and base. They found that f could then take on a value 
of 3.5 and concluded that the decay of the solar flares in- 
volved a large downward velocity at temperatures between 
5 x 10 5 and 10 7 K. However, for the flare loop to drain at 
a reasonable rate (i.e., in a few minutes), they invoked an 
order-of-magnitude change in the flux tube width from top 
to bottom. In support of the idea of flare cooling via super- 
sonic downflows, we note that Dere and Cook (1979) have 
observed large pressure differences between the top and base 
of a flare loop. However, BCS has not observed obvious red 
shifts of the order of magnitude required during flare de- 
cays. This could be due to a number of factors: the down- 
flowing plasma could be at a lower temperature than BCS 
can see; and further, it is not an absolute spectrometer and 
so all the line shifts are relative. 

One common feature of the above studies is that large 
changes in the cross-sectional area of the flux tubes with al- 
titude are invoked. The magnetic field is concentrated in the 
lower atmosphere to strengths of 1-2 kG (Stenflo 1976). We 
have only a rough idea of the actual coronal magnetic field 
strength. Estimates are typically in the range of a few hundred 
gauss (Schmahl et al ., 1982), giving a magnetic divergence 


of about 3 or 4. It also seems that most of the divergence 
occurs in the lower atmosphere (Stenflo 1976), i.e., in a 
region in which the physical validity of the above models 
is doubtful. It is therefore difficult to accept that flux diver- 
gence plays a significant role in the flare cooling process, 
and so the supersonic cooling model does not seem viable. 

Finally, the decay phase of flares has been studied nu- 
merically (Antiochos and Krall, 1979; Antiochos, 1980; 
Nagai, 1980; Doschek et al . , 1982; Pallavicini etal ., 1983) 
with similar results. The first four of these papers all ob- 
tained large scale, supersonic downflows in the decay phase. 
The velocity of these flares was typically 50 - 70 km s _1 , 
implying that for them to be supersonic, the temperature had 
to be less than about 2 x 10 5 K. It is difficult to see how 
such flares can explain the slope of the emission curve up 
to temperatures of 10 7 K. The Pallavicini simulation, al- 
though run for long enough, does not discuss the dynamics 
of such flows. It therefore appears that the physical mech- 
anism responsible for the T 3 5 slope of the DEM is still un- 
clear. Supersonic downflows seem to come closest to 
explaining this phenomenon, but the problem of the ex- 
tremely short draining timescales must be resolved. One 
potential solution is to relax the fluid approximation, as has 
been done in quiet-Sun studies by Shoub (1983). 

Another problem in studying flare cooling is that of fill- 
ing factors. As has been stated in Section 5.3.2, the possi- 
bility of filling factors of the order of 1 % or less means that 
the electron density can be wrong by an order of magnitude 
if derived directly from the emission measure and estimated 
volume. The obvious immediate effect of this is to reduce 
the radiative cooling time and increase the conductive cool- 
ing time by the same factor. This effect can easily be allowed 
for in the cooling time formulae. However, such an adjust- 
ment may be too simplistic. It assumes that any fine struc- 
ture giving rise to such a small filling factor is composed 
of long, thin strands which still communicate to the lower 
atmosphere and still lose energy by conduction, although at 
a reduced rate. An alternative picture is that the filamenta- 
tion gives rise to small, isolated regions in the corona that 
do not communicate with the photosphere and so are unable 
to lose their energy by conduction. An example of this is 
the magnetic islands that arise in the tearing mode (Van 
Hoven, 1981). These small parcels of plasma can then only 
cool by radiation and so give cooling times inconsistent with 
the values derived from the equations simply adjusted for 
the filling factor. Clearly, our understanding of the coronal 
fine structure during flares will be important in determining 
the bulk properties of the decay phase of flares. It is worth 
noting that the three SMM flares (April 8, April 30, and 
November 5) in which we can compare the density derived 
from X-ray line ratios with that from the emission-measure 
volume approach are all large volume flares. An interesting 
goal for the renewed SMM operations should be to make 
equivalent observations of a compact flare. 
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A major source of error is that the behavior of a plasma 
in a temperature gradient is not well understood. The form 
of the conductive cooling time derived earlier relies on the 
Spitzer-Harm (1953, hereafter SH) description of heat con- 
duction, which descibes a plasma only slightly displaced from 
a Maxwellian. In fact, under not very rigorous circumstances, 
the SH form of the heat conduction coefficient fails, and in 
flare plasmas it may fail badly. Rosner et al. (1984) have 
given a concise review of these failings, and we now 
paraphrase their arguments as applied to flare plasmas. The 
SH description of heat conduction is valid only if the mean 
free path (X) is much smaller than the temperature scale height 
[L(T)], where the ratio of the two, R, is defined as 

R = \/L(T) = 10 10 F c /(N e T 3/2 ) < <1. (5.3.17) 

The inequality is strong in the sense that if R = 1/30, 
the SH calculation overestimates the heat flux by a factor 
of 2 (Bell et al ., 1981). In the case where R = 1 (the free- 
streaming limit), the error is at least a factor 10, and labora- 
tory results suggest that the error may be even larger. For 
a flare heat flux of 10 10 ergs cm -2 s' 1 and a density of 10 12 
cm 3 at a temperature of 10 6 K, i.e. numbers characteristic 
of the lower corona, we find R = 0.1. Care must be taken 
in numerical work to check how well the SH method models 
the heat flux; Craig and Davys (1984) outline a method that 
involves a gradual reduction of the heat flux as one moves 
away from the SH regime rather than the sudden jump used 
by many other authors. 

Let us quantify these ideas by considering the five prime 
flares. Table 5.3.5 shows the principal parameters of each 
flare. The cooling times due to conduction and radiation with 
a filling factor of unity have, in some cases, large uncer- 
tainties as a result of lack of accurate information about the 
loop length. The other timescales given in Table 5.3.5 are 
as follows: r mot i on , the time constant of the temperature de- 
cay due to enthalpy flux as proposed by Veck et al. (1984) 
and defined by Equation 5.3.9; 7d ra j n (= L/C s ), where C s 
is the sound speed, is the time taken for the material to drain 
from a loop assuming a large pressure gradient and super- 
sonic flows. 

First, note that r^ m is much shorter than both r 0 b s and 
^motion * The differences between these timescales is that the 
former allows for supersonic motions and the latter do not. 
Thus, the fact that Tdrain < < T obs suggests that the super- 
sonic downflows postulated by Antiochos and Sturrock 
(1982) are not present. Consider the loop to be filled with 
a plasma for which the filling factor is unity and the thermal 
conduction is described by SH. Then the May 21 and Novem- 
ber 5 flares both require long-term energy deposition, since 
the difference between their calculated and observed cool- 
ing times are so significant. The other flares, however, would 
not appear to require further energy deposition. This is con- 
sistent with the predictions of Moore et al. (1980). 


Table 5.3.5 The Principal Parameters and Timescales 
of the Prime Flares 


Date 1980 


Apr 8 May 21 

Jun 29 

Aug 31 

Aug 31 Nov 5 

Time (UT) 


03:07 

21:00 

18:21 

12:49 

12:52 

22:33 

T e (10 6 K) 


8.9 

10.0 

20.0 

12.9 

14.8 

8.9 

11 ,. (10"cm 3 ) 

0.4 

1.0 

3.5 

3.2 

3.8 

1.4 

L (10 9 cm) 


1.6 

3.2 

6.3 

0.5 

0.5 

1.6 

T c 

(s) 

250 

950 

2400 

50 

46 

1280 

T r 

(S) 

4500 

2000 

1150 

800 

750 

230 

T motion 

(S) 

4600 

3960 

1163 

720 

600 


T drain 

(S) 

67 

95 

20 - 80 

14 

13 

52 

T obs 

(S) 

1200 

12000 

1300 

25 

45 

600 


If the filling factor is now set to 0.01 but still with SH, 
the radiative cooling times decrease by an order of magni- 
tude. Then the May 21 and November 5 flares cool more 
quickly than before, but the June 29 flare would also require 
long-term heating, since it cools very effectively by radia- 
tion. Neither of the August 31 flares would require any heat- 
ing, but their cooling is now dominated by radiation rather 
than by conduction. This is a function of our choice of fill- 
ing factor, since a filling factor of 0. 1 to 0.0001 would make 
this conclusion more doubtful. However, if thermal conduc- 
tion is inhibited but the filling factor is 1.0, then the August 
31 flares may need long-term heating, whereas the June 29 
flare does not. 

Even with these simple approximations, the overall pic- 
ture of flare decay is confusing. All that can be concluded 
is that if the observed decay time is longer than the radiative 
cooling time for a filling factor of unity, then long-term 
energy deposition is required, since any smaller filling fac- 
tor will simply make the plasma cool more quickly. 
However, if the flare appears to cool principally by conduc- 
tion with the right timescale, care must be used since other 
factors can invalidate this conclusion. One must also be care- 
ful when invoking long-term heating, since the need for this 
can be removed by a combination of physical effects. Peres 
et al. (1984) found that the computed parameters decayed 
too rapidly when their model was compared with the obser- 
vations, and so concluded that long-term heating was re- 
quired. They found a very simple function that reproduced 
the observed light curves quite well. However, it is equally 
possible that the extended decay could be caused by a filling 
factor effect or by a conductive heat flux limitation. We can 
therefore only conclude that the question of how a flare 
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plasma cools is less clear than it was before. In reality, our 
knowledge of the microphysics of plasmas is so poor that 
we must consider a large number of possible models before 
making any positive conclusions. 

5.3.6 Summary of Conclusions 

Let us address the five conclusions of Moore etal. (1980) 
and suggest how future theoretical efforts could clarify them 
further. 

(1) We find that, under the assumption of filling factors 
equal to one, the densities obtained in these five flares are 
consistent with those of Moore et al. However, the SMM 
results suggest that the particles are contained within fila- 
mentary structures occupying no more than 1 % of the ob- 
served volume (de Jager et al . , 1983; Wolfson et al . , 1983). 
Thus the densities quoted from the Skylab results are too 
small by at least an order of magnitude. The large uncer- 
tainties in the density and possibility of fine structure dic- 
tate our response to the remaining four questions. 

(2) Depending upon the flare parameters, either conduc- 
tion or radiation can dominate the cooling process. Mass mo- 
tions may also play a role in energy transport, a point not 
considered by Moore et al. We therefore disagree with their 
conclusion that conduction and radiation are generally equally 
important. For small filling factors, radiative cooling will 
dominate. 

(3) In the May 21 flare, continued heating is needed, con- 
firming the conclusions of Moore et al. The understanding 
of how this heating can occur has improved on the basis of 
the work by Forbes and Priest (1982, 1983a, b) and Cargill 
and Priest (1982, 1983), which is discussed in Chapter 1 of 
this report. 

(4) The question of whether compact flares need long- 
term heating is now open again as a result of the filling fac- 
tor problem. We do not know what type of filamentary struc- 
ture exists in such flares. One possible way of heating such 
flares is via turbulence (Bornmann, 1985c). 

(5) The chromospheric evaporation scenario proposed by 
Moore et al. seems to be confirmed by our work on these 
five flares. 

Thus, two of the original conclusions of Moore etal. have 
been unambiguously confirmed by our work. Unfortunately, 
the SMM data have confused the other three; but they still 
could be true. A crucial issue for theorists and observers is 
to determine the nature of the fine structure in flares. This 
would appear to hold the key to our understanding of the 
decay phase. ^ 
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5.4. RELATIONSHIPS AMONG THE 
PHASES 

C.J. Crannell and H.S. Hudson 


5.4.1 Introduction 

The overall flare process involves phenomena we have 
characterized as the “impulsive” and “gradual” phases, fol- 
lowing the X-ray signature first recognized by Kane (1969). 
In addition, evidence exists for a pre-flare phase in some 
flares, and recent SMM data have shown that a post-flare 
phase, in which extensive and energetically important coronal 
activity occurs, may also exist. The data to describe the pre- 
flare and post-flare phases are insufficient to place them 
properly into an overall picture of the energetics, aside from 
noting that these phases may indeed be significant from the 
energetics point of view. In this section, therefore, we review 
what is presently known and comment about the possible in- 
teractions among the flare structures involved. 

5.4.2 Relationship Between Impulsive and 
Gradual Phases 

The distinction between the impulsive and gradual phases 
of a flare was originally made by Kane (1969). The ener- 
getics relationship between these phases has been contro- 
versial from the beginning, when Kane and Donnelly (1971) 
showed that the large energy in 10-1030 A bursts correlated 
well with the energy inferred for 10 to 100 keV electrons, 
assuming non-thermal bremsstrahlung as an explanation of 
the hard X-ray bursts. This was the first real evidence that 
particle acceleration during the impulsive phase could have 
energetically significant consequences in the chromosphere, 
the source of the EUV flashes. We now have far better data 
with which to examine this question quantitatively, and this 
subsection deals with the investigation of the energetic rela- 
tionship between the impulsive and gradual phases. 

That the impulsive and gradual emissions are related can 
be seen in Figures 5.4.1 and 5.4.2, in which the peak count- 
ing rates of hard X-ray bursts (Figure 5.4.1) and the total 
hard X-ray counts (Figure 5.4.2) are plotted versus the peak 
Ca XIX soft X-ray counting rate. These scatter plots include 
all events observed with HXRBS and BCS which have both 
a peak hard X-ray counting rate greater than 100 counts s _1 
and a peak soft X-ray counting rate greater than 40 counts 
s~L The total hard X-ray counts (TOTAL) are better cor- 
related with the peak Ca XIX counting rate (BCS) than is 
the peak hard X-ray counting rate. This result supports the 
conclusion of Neupert (1968) that the gradual soft X-ray 
emission resembles an integral of the impulsive hard X-ray 
emission. 

The significance of these results must be evaluated in view 
of the Big Rare Syndrome (BFS) identified by Kahler (1982). 
He found, quite simply, that bigger flares are bigger at all 
wavelengths. Quantitatively, the BFS is manifest as corre- 
lation coefficients of approximately 0.48, with a range of 
0.3 to 0.65, between widely diverse parameters. The cor- 
relations shown in Figures 5.4.1 and 5.4.2 clearly indicate 
a closer relationship than would be expected from the BFS. 
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